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c 

ABSTRACT , This report discloses the full FORTRAN IV pro- 
gram (for an IBM 7094 digital computer) for the automatic 
design of arbitrary TEM quadrature couplers and differential 
phase shifters. The design is completely general in the 
sense that any number of coupled quarter wavelength sections, 
any nominal coupling value (or phase shift), and any design 
bandwidth of operation can be realized, although the user 
may specify any maximum coupling value in the design. This 
last degree of freedom in the specification to the machine, 
which is all-important in the physical realization of a 
theoretical design, is met by finding the proper number and 
types of tandem coupled junctions in the solution. This 
computer solution will be optimum. It will find the least 
ripple for the required bandwidth of operation for any 
(input) design complexity. 

This report shows a complete flow chart of the total 
program. It also gives the special subroutines developed 
for automatic plotting of the coupler (or phase shifter) 
frequency-response functions. Various sample input data 
and machine outputs are also included. 
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Section 1. INTRODUCTION 


During the last decade, before the directional couplers 
attained their present peak of popularity, there has truel^ 
been a host of papers on this subject. (Ref. 1-10.) However 
aside from the question of which published designs could, in 
fact, be built in practice, there still existed the need for 
high-quality (i.e., low ripple) extremely broad band couplers 
Furthermore, because of the recently developed antenna feed 
matrices and receiver techniques, the quadrature relationship 
in some couplers and differential phase shifters are now 
applied in a multitude of microwave system designs. 

This paper is primarily directed to the system design 
engineers who require arbitrary bandwidth, arbit.rary ripple 
(or quality) performance, arbitrary but specifically control¬ 
lable construction constraints on devices and still have an 
optimum solution. The paper is prepared for the engineer 
who does not care about the details of coupler synthesis, 
but only in what the solution and performance is for some 
given requirement and some physical constraint. It is beyond 
tb<=> scope of this report to delve into the synthesis of such 
arbitrary couplers and phase shifters.^ 

The explanation given on the comment cards in the begin 
ning of the program were judged to be sufficiently complete 
so that, any reader may be able to apply the program. This 
explanation became necessary when, in order to make the pro¬ 
gram more flexible, it was modified to accept several types 
of input data. With the enclosed examples, there should be 
no difficulty in using the program. 


^Such is the purpose of a future paper: Shelton, J.P. 
and J.A. Mosko. "Synthesis and Design of Wide-Band Equal- 
Ripple TEM Direction Couplers and Fixed Phase Shifters," 
IEEE TRANS, Vol. MTT-14, No. 10 (October 1966). 
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Section 2. PROGRAM LISTING 


The following text is a listing of the FORTRAN IV pro¬ 
gram (main program and subroutines) to design TEM couplers 
and phase shifters. 

A casual glance at the text will reveal that not all 
subroutines are disclosed in detail. Those missing are 
non-essential in the sense that they are used to plot and 
label responses and write comments using the SC4020 plotting 
equipment. If this plotting system is not available, 

SCOUT V, TABL IV, AI CRT 3, etc.., and SUB ROUTINE PLOT would 
not be reproduced. 

Note that tapes 5, 6, and 7 are for input data, written 
output, and punched card output. Also, tapes IS and 18 are 
for written output (comments, labels, etc.) using the 
SC4020 plotted output. 

Finally, the first five data cards, which are shown to 
be read in, are for comments and labeling of plots: these 
are non-essential and therefore can be blank cards or the 
read and write instructions can be removed altogether. 


2 
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CALL PRINTV (12 .DATE , 
CALL LA6LVI2PC.950.25 
RE TURN 
End 
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SUBROUTINE TABIC(CT,HUT,UMAXT , R T ,NT.UOT.»WT,HIPT,TOLT,RlT,DAT| 
C ft«Utt,*CKO) 


THIS RROCAAA MEfARCt TAILtS 07 OESICN DATA CAOURED IT AT.NT.6WT. 
THIS MOIRA* ACSUIRES THt IAX SUIAOUTlNr 

DIMENSION TOL 110) .A! A ISOS ,C (10,50,10) .III (SO) ,AU10> . C T (10,30) , 

2 NUdO.lO) ,NU*(3Q) ,NUU(30) ,DATE(2) .BAND (20) 

10 fOANAT( 3I2H2(,U ,1M,I2,)W) .fAlCP 9.5) 
to rOAMATt SItHK.II .IH.U.IH) ,2X10710.5) 

•II FOftttAT UltHN* .It.rSXtOHCOU^lINl: , 7*.2, 3H 08) 

90S FORMAT (1X2HN:,12,7SX12HRHA3C SHirTs.rs.l.AH DEC) 

si format (txtOHBANO width , iorio.5 ) 

91 FORMAT (1X11HBAN0WIDTH ,10r 9.3 ) 

9* rOAHAT(IXllHAIAAlt, OEC.lOr 9.5) 

94 rOAMATI1XI1HTOLERANCE 4.10F 9.5) 

99 rOAMAT(1X19M4 READS IN DECREES) 

95 FORMAT (1XIIMACF. LENGTH. lor 9.5) 

97 rOANAT (1X14M* REAOS IN DB) 

99 rOANAT(IXIO hAIRRLC, DB.IXIOP 9.5) 

100 rOANAT <1H ) 

101 rOANAT C1 Hi) 

102 FORMAT(3X7HDATE * .2X2AS) 

AADs 57.295 7 795 

CALL FR9UTV(5) 

If(NCNO.CT.O) CO TO 24 
If(KRUN) 350.350,24 
24 CONTINUE 

BAND(KAON ♦lls BWT 
TOL(KAON)=TOLT 
AIR(KHUN)sAIR T 
RL(KRUN)salt 
OO 45 isl'HNAXT 
ITsNUT(I) 

NU(S,ARUM)sIT 

0 O 45 IAs 1 ,1T 
C(I *IA.KAUN)sCT(|»IA) 

45 CONTINrl 

IT(NENO.CT.O) CO TO 200 
If (NMAX-NMAXTI 47,47.49 
47 AMAXsMNAXT 

49 If(KAUN-10)50,200,200 

50 ir ckt.ne.kcasC) (O to zoo 

If (NT.NE.NTCST) CO TO 200 
TlsrLOAT(NT)400 T 
IfCT1-TCST9.9S) 200,SO,55 
55 If(T1-TCST41.02) SO,SO,200 
SO At TUAN 
200 CONTINUE 

WAITE(IS,101) 
iriftCASC) 210.210,250 
210 CONTINUE 

WAITE(IS.100) 


l 




350 TCSTxfLOAT (NT) *i 
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Section 3. FLOW CHART OF THE PROGRAM 


Using an automatic program with the SC4020, a lineal 
flow chart was prepared of the previous program. Because 
of the lineal display and the number of pages that are 
required, wherever possible, two "frames'' were placed on a 
single page in the report. Thus, on page 19, for instance, 
the heavy broken line represents the proper sequence of 
events. In the absence of this broken line, it will be 
understood that flow lines at the top on one page directly 
connect to lines at the bottom of the immediately preceding 
page. 
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confute 17CMAX 


411 

WftlTC(«,•»>I 
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Section 4. SAMPLE APPLICATIONS 


Example A 
Type 1 Data 

If it is desired to design a 17:1 bandwidth, -3 db~ 
coupler, which has a normalized Eoe ± 1.8343 limitation on 
it, one needs the Chebyshev antenna distributions for 15 
elements at 19.6 db side lobes: 

1.0, 0.962, 0.890, 0.789, 0.663, 0.536, 0.837 

Assuming that three iterations are allowed (1% toler¬ 
ance (.01) is desired) and that neither development (TDEV=0) 
nor printed performance (X?RINT=0) data are required, then 
the punched input data cards are as shown in Fig. 1. The 
total output data is shown in Fig. 2, 3, and 4. 

The computer execution time was 5.93 seconds. 
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|0 0 010 0 010 0 0)0 0 0|0 0 0|0 0 0|0 0 010 0 01 0 010 0 010 0 010 0 010 0 010 0 CIO 0 010 0 010 0 010 0 010 0 010 0 010 0 010 0 01!) Q 010 0 010 0 0 0 0 0 0 0 

a i.oooi ic. 


IQ 0 Q!0 G 010 0 Gio 0 010 4 V 0 OiO 0 010 4 OlO 0 OlO 0 010 0 $10 0 0!0 0 CIO 0 010 0 010 0 010 0 010 0 010 0 G|0 0 9|0 0 013 0 010 0 0|0 0 0|Q 0 0 0 0 0 0 9 

1 ? -n 3 t.r.COOO 45.00000 O.OIOOO 1.8;.*3*+ 17,CO000 


00000fl a 000n0000 4i,j,, 0a000 ,4,44 10 3a 1 2 3 4 5 0 , j0 4 ' ! 4 000 0-0 000l50!n!000 4 4 j j 44 00900000 0 je 0030a00000 

i 2 > « i I i » » u ii n o M is '» n u il Miiu;);iiliHi aaii|!i]i ii n nx bum hi.h i nuts iim 11 n ■; u i> n it » inn 

T ’ m 1 rmTTr r i i minimi TTfiT 11 TTrTili i mm FIT muimii f TTTh i m m m 


222222222 ? 22 ? 22222222222222222222222222222 ? 2222222222222222222222222222?22222222 
3 3 3 3 3 3 3 3 3 » 3 3 3 3 6 3 3 3 3 3 3 3 3 3 6 3 3 3 3 3 3 3 3 3 • 3 3 3 i 3 3 3 3 3 6 3 * 3 " 3 3 3 3 3 * i 3 3 3 3 3 3 3 Z 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 

4444444444444444444444 '444444444444444444444444 '* *444444M4<444444444<44<44444 M 

5 5 5 5 5 5 5 S 3 5 5 5 5 5 5 5 5 4 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 * s 5 5 5 5 5 5 5 5 b 5 5 5 5 S 5 5 5 5 

1 

6 6 6 666668666666 $ 6666666666 B 6666 G 5 6 6 6 6 6 6 6 6 6 6 6 6 6 5 6 S 6 6 6 6 6 6 G 6 S 6 6 6 8 S 6 & 866 S 6 6 6 G 6 B 6 BBS 6 

7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 K7 7 7 7 7 7 7 7 *77 t 7 7 7 7 7 7 7 7 7 7 7 7 7 ? 7 7 


8888S3S & 8.8 8 8 8 8 8888SS8B8 8 8 5 8 8 8 3 8 8 8 3 8 s|s 8 8 8 S 8 8 8 8 818 8 8 S 8 8 8 8 t 8 8 » 8 8-8 8 8 8 8 6 8 8 8 8 8 8 8 8 

ii ii 


9 9 9 9 9 ? 8 9 9 9 8 2 9 9 9 9 3 9 9 9 8 5 9 9 9 9 9 9 9 9 9 9 9)9 9 9 9 9 9)9 9 9|$ $ 9 9* 2)3 8 9 9 8 919 9 9 9 9 9 9 8 9 3 9 9 9 S 9 9 0 9i9 5 9 9 9 S 9 9 

i 2 jl« i «n * • iio it ufu w h ii ’i umiiIjk) jt!;i» jiwmj: n j;»»<»» ins <n.c. >^stsi wv'.'l'iKstwnttwiMi’WMK’sifHs if •••■•>< * " tto 


FIG. 1. Input Data Cards 


Assuming this to be the. first design problem for the 
machine, then these input data would be preceded by five 
cards: 

1) Instable B f s 

2) Failed 

3) 

4) Theta (degrees) 

5) Coupling (db) 

with the reference 


12345 


Note - Card 3 was left blank since the plots (see Fig. 6, 7 
and 8) were not desired to be titled on the top. 
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17.000 TO 1.0 BANDWIDTH COUPLER 


MS4I4 

9.9417ft 

1.99494 

0.9417ft 

1.94797 
0.49179 

1.94001 

0.9119ft 

1 .39440 
0.20798 

t.14580 
0.19479 

1.04944 
0*04949 

1••9414 
C.9417ft 

£.>2990 
O.1962ft 






1.99494 

ft.9417ft 







1*9194? 
ft.Oltftft 








FIG. 2. Response Curve and Design of Coupler 
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FIG. 3. "Written” Machine Output 
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FIG. 3. (Continued) 
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Example B 
Type 2 Data 


To improve further on a coupler design, one on which 
data are available but in terms of reflection coefficients 
Isuch data are part of the outpr from previous calculation ’ 
e,g., Fig. 4), assume that the uata are as in Fig. 5 


I’ = 0.4989, r 2 = °- 1536 > r 3 = °- 0881 > r 4 " °- 0575 > 
r 5 = 0.0403, r 6 = 0.0277, r 7 = 0.0303, 


with six allowable 
and 1% tolerance, 
Fig. 6. 


iterations, required development data, 
the input data will then be as shown in 


The output data for this computation are shown in 
Fig. 7—11. 


The computer execution time was 6.30 seconds. 










NAVWEPS REPORT 9048 


1 » » ™ '1 »* » w P 2 l_ *11™ J?_*: “J 3 74 r>W7/ Tin 30IJI 32 33 34 3 3 » 3/ 3131 40 «» 4i «3 54 IS 4( 47 4! 4* SO St SJ S3 S4JS5 Sf S7 S» 5160 68 6? 6J|64 6S It 6J M 69 70 .1 72j?J *1 7S 76 77 7) 79 10 

innnrm TTij? fill 11 mm i \ 11 iTiimrnmn uTffi n 11:11 iTil mfi nil m n n 


2 2 2 2 7 2 J 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 7 2 7 . 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 

3 3 3 3 3 3 3 3 3 0 3 3 3 3 3 3 3 3 3 3 * 3 3 3 3 3 3 3 *3 3 3 3 3 3 3 3 3 3 “ 3 ” 3 “ 3 3 3 3 3 3 " 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 
44 4j4 41444444444444 *4 4|4 44444444444444444444 '4 '44444444444444444444444444444444 


-f ------j---- j - 

44 4j4 44444444444444' 44444*44444444444444444'4*44444444444444444444444444444 

5555555553555555555 555555555555555555555555555555555555555555555555055555 
6 E E 6 5 E £ S S S 8 £ 6 E S 6 6 6 S 6 B 6 6 5 5 6 6 G 5 5 6 S 8 6 5 6 6 6 6 5 6 5 6 a 6 6 5 5 6 6 6 S 6 6 6 6 6 fi 6 6 6 6 6 E 6 S 6 6 Rtf 6 6 6 6 


5 5 5 5 5 5 


7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 717 7 7 7 7 7 7 7 7 7 7 / 7 7 7j7 7 / 7 1 7 7 I , 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7|7 7 7 7 7 7 7 7 7 7 7 

58688S83E 8 8 8 6 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 sj* 8 8 8 8 8 8 8 8 1 8 1 8 6 8 8 k 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 S 8 8 8 8 

9 5 9 8 8 3 9 8 9 8 3 3!9 8 9 3 9 9 9 3 9 9 $ 318 9 9 9 9 9 5 9 3l9 3 3 9 9 9 9 S 9 9 9 9 3 9 9 8 S S 9 9 8 9 9 8 9 9 9 9 9 9 9 9 9 8 9 SiS 9 5 9 9 9 9 9 

4 s i i i i ;'8 «nj :• rt|u i» ti if 2c 2'..;: < n ?s ;>! i•* j< j\ jilr # s ; v ' 4 M 44 4; sc v s; sj si ss« sj si s«« n «*» v« ss $6 * 6 ? wio n * 1*5 s a •*? .s 910 

J#<i^2!4* ' 

FIG. 5. Inpuc Data Cards. 
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FIG. 7. Response t rve of First Tandem Coupler of Solution 
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THETA (DEGREES) 

FIG. 8. Response Curve of Second Tandem Coupler of Solution. 
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FIG. 9. Third Tandem Coupler's Response. 
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Example C 

Type 3 Data (Phase Shifter) 


Let it be required to check and improve, if necessary, 
a phase shifter design that is in terms of even mode im¬ 
pedances: specifically, the 17:1 bandwidth, 90 degree phase 
shifter of Ref. 12. Thus, the input data are as in Fig. 12 
with the output data as shown in Fig. 13-15. 

The computer execution time was 2.85 seconds. 


Example D 

Type 3 Data (Coupler) 


Entertain the request for a 7.50:1 bandwidth, -3 db, 
re-entrant nine-section coupler with a normalized Zoe = 4.5 
for the maximum coupled section. Note Ref. 10, which con¬ 
tains such a coupler design, except that Zoe = 4.93133 and 
that the bandwidth is 7.365:1. Therefore, let this published 
data be the starting point of the desired data. Thus, the 
input data cards are prepared according to Fig. 16 with the 
output data shown in Fig. 17-19. 

The computer execution time was 7.65 seconds. 
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17.000 TO 1.0 BANDWIDTH PHASE SHIFTER 

1.0*414 1.49434 1.41434 S.ftS*** 1.S7341 1.1*034 
0.4417* ft.*411* 0.1417* 0.4*130 9.30701 0.17*31 

1. *3434 1.43.3* l.t4**Z 

0.1477* 0.1417* !>.*»»•• 

1.414.14 1.1*0*0 

0.*417* 0.1*4(1 

I.SI40* 

0.414*7 


FIG. 13. Response Curve and Design of Phase Shifter. 
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FIG. 14. ’’Written” Machine Output 
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FIG. 18. (Continued) 
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FIG. 18. (Continued) 
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Example E 
Type 4 Data 


Should it be desired to calculate a coupler similar to 
the presently completed one of example D except for Zoe =3.0, 
then one could take the output data of example D (Fig. 19) 
to be the input now. However, one can save the punching of 
T's, Zoe's, etc., by punching CONT = 5. on the first data 
card. This simply saves the design data from the previous 
design to be the new starting data. Thus, Fig. 20 was pre¬ 
pared as input and Fig. 21-23 are its output. 

One should use the Type 4 data input with some discretion. 
One would not, for instance, start the design of a 5-section 
phase shifter based on a previously completed design of a 
6-section phase shifter.3 However, one can start a solution 
by using a presently completed solution of basically the 
same device except changes in Z max. and BW, or both. 

Obviously, the greater the differences between the known 
and the desired designs, the longer the iteration process 
and the less efficient the method of design. 

The computer execution time was 3.72 seconds. 


interesting to note that in preparing the input 
data tor example C above, an error in copying data was made 
resulting m a deletion of a coupled section of Zoe = 1.8343 

c^k 0 the -! 1 u C tu ne J Star l ed the desi § n of a 5-section phase 
shifter with the data from basically a 6-section coupler 

(see Fig. 24 for actual input data cards). To see how the 
machine performed under this input, see Fig. 25 where the 
performance response s maxima/minima (peaks/valleys) are 
^ (between 0=0 degree to 0 = 90 degrees), also see 
Fig. 26 for final response plot. Thus, the initial response 
does not have the proper" number of maxima/minima, but 
after one iteration, this has been rectified. No complete 
analysis has been made to specify, in general, how poor a 
starting solution can be and still have the machine converge 
on a meaningful result. 6 
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FIG. 25. End 
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FIG. 26. Final Response Curve and Design Data of Phase Shifter. 
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